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Abstract. This work studies some aspects of the computational power 
of fully asynchronous cellular automata (AC A). We deal with some no- 
tions of simulation between ACA and Turing Machines. In particular, we 
characterize the updating sequences specifying which are "universal", i.e., 
allowing a (specific family of) ACA to simulate any TM on any input. 
We also consider the computational cost of such simulations. 
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1 Introduction 

Cellular Automata (CA) are a computational model widely used in many sci- 
entific fields. A CA consists of identical finite automata arranged over a regular 
lattice. Each automaton updates its state on the basis of its own state and the 
one of its neighbors. All updates are synchronous. 

CA are particularly successful for modelling real systems [5]. However, many 
natural systems have a clear asynchronous behavior (think of biological processes 
for example [H]). Asynchronous CA (ACA) have been introduced in order to be 
able to more closely simulate these systems. Roughly speaking, they are CA in 
which the constraint of synchronicity has been relaxed. 

According to which updating policy is chosen, the behavior of the ACA un- 
der consideration can be very different. In literature, several policies have been 
considered (purely asynchronous |15| . a- asynchronous |18| . etc.) 

In this paper we consider a fully asynchronous behavior in which, as in a 
continuous time process, two cells are never updated simultaneously (or, equiv- 
alently, one and only one cell is updated at each time step) |6I19) . Of course, 
the evolutions of such ACAs depend on the sequence of cells that are updated 
at each time step. 
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It is well-known that CA are capable of universal computation (see [T7] for 
an up-to-date survey). Various mechanisms have been introduced to show how 
AC A can emulate classical CA or circuits |15I16I11I21| . As a consequence of 
these results, Turing universality of AC A is also proved. 

In this paper we focus on the way by which an ACA simulates a TM. We 
analyse two different modes: strict simulation and scattered strict simulation. 
Roughly speaking, strict simulations pretend that the ACA exactly reproduces 
the steps of the TM (up to some encoding) admitting at most that some time is 
wasted between two steps of the TM. The second mode is essentially the same as 
strict simulation but considers the case that only a subset of cells can be used to 
perform the simulation, the others being "inactive". We characterize all updating 
schemes allowing these two simulation modes. Moreover, we show that in both 
cases, the time slowdown due to asynchronicity is quadratic w.r.t. the running 
time of the TM under simulation. 



2 Basic Notions 

Let A be an alphabet. A CA configuration (or simply configuration) is a fimction 
from Z to A. A cellular automaton is a structure (A, r. A) where A is the alphabet^ 
r is the radius, and A : A^^^^ — > A is the local rule. The local rule is applied 
synchronously to all positions of a configuration. In other words, the local rule 
induces a global rule / : defined as follows: 

Vc G A^,\fi G Z, f{c)i = A(ci_r, . . . ,Ci+r) 

Where A^ is the set of all configurations and for each i G Z, the image c{i) 
is denoted by Let A be a local rule of radius r. Consider now the following 
asynchronous updating of a configuration. At each time i, / is applied on one and 
only one position. We denote by 6 the sequence {Ot)t>a whose generic element 
9t is the position which is updated at time t. 

Definition 1. A fully asynchronous cellular automaton (ACA) is a quadruple 
{A, A, r, 6) where A is a finite alphabet, A : A'^^^^ ^ A is the local rule of radius 
r G N and = {0t)t>a, with 6t G Z is a sequence of cell positions. 

Every ACA A — {A, A, r, 6) induces a dynamical behavior described as fol- 
lows. The evolution of any configuration c G A^ is the sequence of configurations 

{f{c)j\c),f\c),...j\c)j'+\c),...}, 

where /"(c) — c and, for any t G N, /*^"'^(c) is defined as 



ViGZ, f+^(c) 




.,/*(c),+,) ifi^Ot 
otherwise 



We stress that dynamical evolutions in asynchronous CA depend on the choice of 
the updating policy. Different updating schemas have been considered for study- 
ing asynchronicity in C A settings |15|6|18| . We deal with the fully asynchronous 



situation in which at each time only one cell is updated according to 6. Note 
that even for some relatively simple rules the behavior under fully asynchronous 
updating is far to be simple, see for instance |12I18| . From now on, for a sake of 
simplicity, by referring to ACA we will mean fully asynchronous CA. 

A Turing Machine is a 7-tuple (Q, Z", J", b, 5, go, F)^ where Q is the set of 
states. S (Z r is the input alphabet, F is the working alphabet and b G -T \ is 
the blank symbol. The map S'.QxFt-^QxUx {L, R} is the transition function., 
where L and R denote the left and right movements of the head, ^ Q is the 
initial state and F C_ Q the set of final states (see [7 , for an introduction on 
this subject). An instantaneous configuration c of is a triple {T,q,p), where 
T £ is the content of the tape, q £ Q is the current state of Ai and p £ Z is 
the position of its head. 

A run of Ai on the initial input x £ S* is the sequence Rt — {(Tt, qt,Pt)}ti£i. 
where (Tq, qo,Po) = ("ba;b", qo, 0) {i.e. the symbol b is repeated infinitely many 
times at the left and at the right of a;)for t = 0, and for any i G N, (Tj+i , qt+i , Pt+i) 
is the instantaneous configuration of Ai at time t + I, where T(_|_i is equal to Tt 
except that in position pt in which the symbol {Tt)p^ is replaced by the symbol s 
with (gt+i,s,X) = 5{qt, {Tt)pJ, where pt+i = pt + 1 if X = i? and pt+i = pt - 1 
a X = L. Remark that if in a run it happens that at some time t £ N the state 
qt £ F, then, for any k > t, Rk — Rt, i-c, in other words the computation halts 
on the instantaneous configuration Rt and the output is the non blank content 
of tape Rt. 

3 Simulation of Turing Machines 

It is well-known that CA are a universal computational model according to 
different notions of universality (for a survey see |17|). The main point to prove 
universality is to simulate a TM. Of course, one can apply similar ideas and 
constructions to prove computational universality or computational capability 
of ACA (see for example |15|16|11] ). In this section we would like to precise the 
computational cost of such simulations. 

The basic idea when simulating a TM using an ACA is to act by "extracting" 
first the information about the state of the TM from the current configuration of 
the ACA, and then to operate the TM transition saving information on the ACA 
configuration again. The way of saving the TM state and the way we extract it 
from the current configuration lead to the two following notions of simulation. 

Notation. For a,b £ N with a < b, denote [a, b] the set of integers between a 
and b (including a and b). Given the finite sets Ai, A2, . . . , An, for any element 
(ai, 02, . . . ,an) £ Ai X A2 X . . . X An, define 7TA;((ai, 02, . . . , Oi . . . ,a„)) = ai. 
These projection maps can be naturally extended to work with configurations, 
indeed given a configuration c £ {Ai x . . . x A„)^, for any i £ [l,n] and j £ Z, 
define 77^^ (c) j = 77^. [cj ) . 

Given a configuration c £ A?' and a function -tj) : "L ~^ c^' is the configura- 
tion defined as cf = c^(i) for all i £ Z. 



Definition 2. Let M = {Q, S, F, 6, S, qo, F) be a TM and A = {A, A, r, 6*) be an 

AC A. A strictly simulates M iff A = F x B for some finite set B and for any 
input X € S* of M, there exists a configuration c € A^ satisfying the following 
conditions: 

1. nr{c) = '^bxb^ and nB{c) = "sus" for some u,s G B; 

2. for any time t e N, there exists e N such that 

nr{f''{c))=Tu 

3. for any pair of times ti,t2 € N, 

ti <t2^ t[ < t'2, 

where t[ = min{fe e N : Bpif'^ic)) = TtJ, i = 1,2. 

In other words, an ACA A strictly simulates a TM A4 if its configurations can 
represent in a direct way the tape of M, possibly using an additional amount 
of information (stored in the alphabet B) and some additional time. Relaxing 
the condition on the representation of the tape, the following weaker notion of 
simulation is obtained. 

Definition 3. Let M = {Q, S , F, bj, qa, F) be a FM. Let A = {A,\,r,e) be 
an ACA. A scattered strictly simulates A4 iff A = F x B for some finite set 
B . and for any input x G S* of M. there exist an injective increasing function 
■)/) : Z n> Z and a configuration c £ A^ satisfying the following conditions: 

1. (a) nsic^) = '^sus'^, for some u,s€B ; 

(b) Vi e Z \ V'(Z), iTs(ci) = q, for some q e B ; 

(c) nr{c^) = '^bxb'^ ; 

2. for any time t €N, there exists t' G N such that 

nr{{f'\c)f) = Tt 

3. for any pair of times ti,t2 € N, 

ii < i2 t[ < 4, 
where t'^ = mm{k G N : i7r((/'=(c))'^) = TtJ, i = 1,2. 

A scattered strict simulation assumes that only a subset of cells participates 
to the simulation and the others are somehow inactive. For this reason, in the 
ACA configurations there can be an offset made for example by bs between the 
symbols of the TM tape content. Note that when the function -ip is ip{i) — i then 
scattered strict simulation and strict simulation coincide. 

According to the above definitions, even if an ACA can simulate a TM on a 
fixed input x, it might not be able to simulate the same TM on a different input 
simply because of an inappropriate updating sequence 6. 



3.1 Construction 1. 



Given & TM M ^ (Q, S, T, b, S, qq, F) build a family of ACA Ae = {A, A, 1, 9) 
such that A^TxQxDxC, where D = {L,R} ,C ^ {0, 1, 2}, and the local 
rule A : A"^ A is defined as follows 



A(u, V, z) 



if M = {au,qu,R, = io-v,qv,mv,2) 
and (5((j„,crt,) = {q,(j,m) 
{a.u,qv,R,'2) if u = {a^,qv,R, l),z = (ct^ , , , 0) 
(cr„,g^,mt,,l) if u ((t„, g„, i?, 2), v = (cr^, m^, 0) 
if z = (cr^,(j^,i,l),u = {ay, q^, my, 2) 
and (5((j2,cr„) = {q,cr,m) 
{cFy,qy,L,2) liv = {ay,qy,L,l),u = g„, m„, 0) 
((T^,g^,m^,l) if z = (cr2,g2,L, 2),w = (cr^ , , m„ , 0) 
v otherwise. 
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Fig. 1. Simulation of the first step of a TM using an ACA built by construction 1 witli 
updating sequence 9 = (0, —1, 0, 1, . . .). The ACA and TM times are denoted by t' and 
t, respectively. The arrow points at the current active cell of the ACA. 



Every cell of the ACA contains the symbol of the corresponding cell on the 
TM tape, the state of the TM, the direction of movement of the TM head, and 
a value ^ S C to control the simulation. At TM time t, the cell i with £, = 1 
is the one where the TM head is positioned at time < — 1, i.e., pt-i = i (with 
p_i = —1). During the ACA evolution, at most one cell in whole configuration 
has ^ = 1. If the updating sequence allows the cell i + 1 (resp., i — 1) to be 
updated and the cell i has m — R (resp., m — L), then the cell i + 1 (resp., i — 1) 
changes its state according to the TM rule and its own value of ^ is set to to 
indicate that the information about the head position has to be moved to this 
cell. To perform it, at subsequent times ACA will set the cell witlr ^ = 1 to 2 
and the cell with f = to 1. An example of this behavior is shown in Figure [T] 

In order to (strictly) simulate a TM on input x = xq . . .Xn-i G ACA 
given by the above construction have to start on the following configuration 



Vi e Z, Ci= { 



{xi,q,m,2) if < i < |a;|. 
(b.go.-R, 1) if i = -1. 
^ (b, q, m, 2) otherwise, 



where q and m are an arbitrarily chosen state and movement (since they will not 
be used in the simulation, their choice can be arbitrary). The last point to precise 
is which updating sequences can be used. Of course, that depends on the TM 
to simulate but there are sequences that can be used in "all occasions", they are 
called universal. An updating sequence 6 is universal iff |{i £ Z, 0j = fc}| = oo 
for every k G Z; practically speaking a sequence is universal if any cell is updated 
infinitely many times. 

Theorem 1. An AC A A = {A,X,1,6) given by Construction 1 is Turing uni- 
versal if and only if 6 is universal. 

Proof Consider a TM = {Q,IJ,r,S,h,qo,F) and an ACA A = {A,f,l,9) 
with universal. For any x G S* input of A4, let c be the initial configuration 
built in construction 1. Let us prove that A strictly simulates A4. 

Let Rt = [Tt, qt,Pt) be the configuration of the Ai at time t. We claim that 
for alH G N there exists t' G N such that the configuration c' = /* (c) of A has 
the following properties 

1. iTr(c') = Tt. 

2. 7Tq(c;^_ J = (p-i = -l) 

3. 77c.(c;j = 1 and i7c(c^) ^ 1, for alH e Z \ {pj. 

4. HniCpJ = Rifpt> Pt-i; L if pt < pt-i- 

We proceed by induction. For t = the claim is true by construction {t' = 
and c' = c). Assume that the claim is true for < > 0, i.e., there exists t' such that 
the configuration c' = /* (c) satisfies the four stated properties. Remark that 
nc (Cpj ) = 1 and hence the only cells that can change their value arc at positions 
Pt + 1 01 Pt — 1, depending on the value of I?£)(CpJ. Assume that nD{CpJ = R 
(the other case is similar). Since 6 is universal, there exists t" > t' such that 
^t" = Pt + 1 and for any other t G N cither Of ^ pt + 1 ot t > t" . According 
to the definition of /, at time t" the cell Pt + ^ will become (a, g,m,0), where 
{a,q,m) = 5(i7Q(CpJ, iTr(Cp^+i)) and c" = /* (c). Moreover, no other cell can 
change its content between time t' + 1 and t" — 1. Therefore Uric") = Tt+i and 
-ffQ(Cp^+i) = qt+i i-e. c" satisfies the first and second properties. 

Again, since 6 is universal, there exists t'" > t" such that Of" = Pt and for 
any other t G N either 9t pt + 1 or t > t'" . According to /, the forth component 
of the cell at position Otm is set to 2. Once more, remark that no changes in the 
configuration of the ACA occur between time t" + 1 and t'" — 1. 

Finally, by the universality of 9, there exists i > t'" such that 9i = pt + 1 
and for any other t G N either 9t ^ pt + I ot: t > t. From the definition of /, 
one deduces that the only possibility is that the fourth component of the cell at 



position + 1 in /*(c) is set to 1. Since no changes in the configuration of the 
AC A occur between time t'" + 1 and t — 1, the claim is proved. 
To prove the inverse imphcation, consider the TM 

M = ({te,gL},{O,l},{O,l,b},b,(5,te,0) 

where 5 is defined as follows 



(9,cr) 






(te, 1) 




(9L,0) 




(5(g,cr) 








(te,0,i?) 


(<ZL,0,L) 


(9L,0,i) 



On any input, M writes a symbol 1 on the cell 0, then it writes Is towards 
the right until a blank symbol is reached. When a blank is reached, it moves left 
writing a symbol until a blank is encountered at this point it starts moving 
right writing Is and so forth. It is clear that the head of M. passes through any 
cell of the tape infinitely many times. Therefore any ACA given by Construction 
1 needs an universal updating sequence to strictly simulate it. □ 

The previous result proves that the class of ACA given by Construction 
1 are computational universal but it seems that requiring an universal updat- 
ing sequence involves a considerable time loss (see the proof of the Theorenfl]) . 
The following proposition shows that there exist (carefully chosen) updating 
sequences such that the time loss is acceptable (quadratic). 

Proposition 1. Given a TM M. that executes in time T{n), there exists an 
ACA A given by Construction 1 that simulates M. in time O {T{n)^^. 

Proof. Let Si with i e N be the finite sequence given by all the integers between 
—i and -\-i with step 2 (e.g., S2 = (—2,0,2) and si — (—1,1)). Consider the 
sequence 9 given by the concatenation of the Si sequences: 

g = so(-l)so siSQSi S2SIS2. . . = 0, -1, 0, -1, 1, 0, -1, 1 , -2, 0, 2, -1, 1, . . . 

Clearly, is universal. Consider A = (A, A, 1,6) where A and A are as in Con- 
struction 1. It is easy to verify that every block SiSi-iSi simulates one step of 
A4. The size of the block increases by a (multiplicative) constant 3 for every i. 
The total length of the simulation is then bounded by 

Tin) Tin) 

l+l^o\+J2 2-k.l + k.-il = 2+5] 3z+2 = - (T{nr + -T{n) + 3] = O {T{nf) . 

Remark 1. The time O {T{nf) is the better asymptotic limit for a sequence. 
Indeed, consider a sequence such that the corresponding ACA simulates in time 
0{g{T{n))) C 0{T{ny) any TM which executes in time T{n). We show that 
0{glT{n))) = 0(r(n)2). Consider all the possible movement of the TM head. 
Its position is at time 0, either 1 or —1, at time 1, and so on. Let us introduce 
the graph G = {V,E) where V = {0,...,T(n)} x {~T{n), . . . ,T{n)} and E = 
{{{t,a),{t + l,b)) 1 t e {0, . . . ,T(n)-l} and la-61 = 1}. Any path starting from 



(0, 0) and ending to {T{n), h) with h e {— r(n), . . . , Tin)} represents a possible 
sequence of a TM head. In order to simulate all the head movements, all the 
nodes of the graph must be visited at least once by these paths. Since the graph 
has 0(T(n)^) vertexes, the time 0{g{T{n))) can not be less than 0{T{nY). 

Remark 2. Actually, for any a TM which executes in time 0{T{n)), there 
are uncountably many ACA given by Construction 1 that simulate A4 in time 
0(T(n)^). An infinite set of them is individuated by the sequences obtained from 
the one illustrated in the Proposition 1 "inserting" in it any integer in one or 
more positions. 

The slowdown in the simulation of TM using ACA given by Construction 1 is 
essentially given by the fact that we want a strict simulation and we must keep 
track (among other things) of the position of the head of the TM. Relaxing 
this last constraint brings to a different notion of simulation and to a different 
construction. 

3.2 Construction 2 

Given a TM 7W = {Q, S, T, b, 5, qo, F) build a family of ACA Ae = {A, A, 1, 6*) 
such that A = r xQ X D xC, where D ^ {L,R} ,C = {1, 2} and A : -> A 
is defined as follows 

if li = {au,qu, RA), nr{v) = cr^, 
nc{z) ^ 1 and 5{qu,(Jv) = {q,a,m) 
if z ^ {a2,q^^,L, l),i7_r(v) = cr„, 
ndu) ^ 1 and 6{qz,ay) = {q',a',m') 
otherwise. 

The initial configuration for starting the simulation is the same as the one given 
for Construction 1. As already remarked before, this family does not strictly 
simulate Ai since it does not keep track of the head position. However, using 
similar techniques as in Theorem [TJ one can prove that for 9 — sqSiS2---, the 
ACA (A, A, 1, ^) simulates M on any input with a total running time 

T(n) T(n) 

Y^\s,\ = J2i^ + l) = - {T{nf + 3T(n) + 2)^0 {T{nf) 

j=0 1=0 

where T(n) is the running time of M on the input x. 

3.3 Construction 3 

Construction 1 and 2 assume that potentially all cells of the ACA cooperate to 
the simulation of the TM. Assume now that only a subset of cells participates to 
the simulation and the others are somehow inactive. In this section, we are going 



A(u, v,z) — < 



(cr, g,TO, 1) 
{a',q\m'A) 

{nr{v),nQ{v),nD{v),2) 



to show that the ACA can still (scattered strictly) simulate any TM whenever 

the updating sequence has some specific properties. 

A set S C Z is syndetic if there exists some finite E C 7, such that UnesiS — 
n) — Z, where (5 — n) = {fc G Z | fc + n G S*}. Syndetic sets liave bounded 
gaps i.e. there exists 17 G N (which depends on S) such that for any h E 1,, 
{ft,, /i + 1, ...,/?, + g} n 5 7^ 0. Given a sequence a = {Qijjgz, the support of a is 
the set supp(q;) = Uigz {ai\. 

Notation. To shorten up the notation in what follows, given an ordered se- 
quence of states it^^''^ . . . ,u^°\ . . . ,u'-^\ denote Enik) = {i G [l,r] | 77c (u^*^) = k} 
and similarly 7i;i(fc) = {i G [— r, —1] | 7Jc(u^'^) = k}. Finally, denote jij = m.mEii{k) 
if Enik) ^ ID and Jl = max EL{k) if EL{k) ^ 0. 

Given a TM 7W = (Q, S, E, S, go, F) build a family of ACA Ae = {A, A, r, 9) 
such that A = E X Q X D X C, D = {L,R}, C = {0,1,2,3}. The local rule 
A : j4 is defined as follows 



i-r) 



,(0) 



,ir)\ - 



(a, q, TO, 0) 



{ay,qy,R, 2) 



if £;l(1) ^0, 
u<^jL) = (cr„,g„,i?, 1), 
u^°^ = {(rv,qv,my,2), 
and (5(g„,cr„) = {q,u,m) 

if = {ay,qy,R, 1), 
i^fl(O) ^ 0, 

and = (fT^,g2,TO2,0) 



if £;l(2) ^0, 
(a-„,g'„,TO„, 1) u^-'^) = (cr„,(7„,7?, 2), 

= (CTt,,qt,,TO„,0) 



(a, q, TO, 0) 



(cr„,g'j;,L,2) 



ifii;fl(i)^0, 

u^in) = (a„g„L,l), 
= {ay,qy,my,2) 
and 5(fe,fT„) = {q,a,m) 
if = (cr^,g„,L, 1), 
i^L(O) 7^ 0, 

u^:>l) = (CT„,g„,TO„,0) 



if (2)^0, 
(cT^,g„,TO^, 1) u^-'^) = (cT^,g^,L,2), 
= (cr„,g„,TO„,0) 



,,(0) 



otherwise. 



In order to be able to (scattered strictly) simulate a TM on input xo ■ ■ ■ Xn-i S 
S*, ACA given by the above construction have to be started on the following 



configuration c 



if i = a_i 
if i € [Q!o,a„-i] 
if i e supp(a) \ [ao, Q;„_i] 
otherwise 

where a is a subsequence of 9 such that ao < ai < . . . < a„, and q Cz Q, m E D 
are arbitrarily chosen. Clearly, the whole construction (and hence the simulation) 
depends on 9 and its subsequence a. The following result characterizes them. 

Theorem 2. An ACA A — (A, A, r, 9) given by Construction 3 scattered strictly 
simulates any TM on any input if and only if 6 contains an universal subsequence 
a such that supp(a) is a syndetic set. 

Proof. For any TM M ^ (Q, U, F, b, S, qo, F) consider an ACA A = {A, A, r, 9) 
given by Construction 3. Assume that A scattered strictly simulates Ai. First 
of all, let us prove that |supp(6')| = cxd. Indeed, if |supp(0)| < oo, only a finite 
number of cells can be used for simulation and therefore, according to Condition 
2 of Definition [3] only a finite portion of the tape can be simulated. Choose 
a subsequence a oi 9 such that ao < ai < . . . < a„ (this is possible since 
|supp(0)| — oo). By contradiction, assume that no a is universal. This means 
that the set of cells that can be updated infinitely many times is finite or empty. 
Without loss of generality assume that it has finite cardinality and j > is the 
maximal of its elements (the case j < is similar). Let k be the index of last 
occurrence of j in a. Consider the TM 7W from the proof of Theorem [T] on the 
empty input. Since, for all t > k, /*(c"))j = c", Condition 2 of Definition [3] is 
violated. 

Now, always by contradiction, assume that there exist universal subsequences 
of 9 but none of them has a syndetic support set. This means that there are larger 
and larger sets [a,b] C N not contained in supp(Q;). Choose one of them such 
that b — a>r and let Ai be the TM which on the empty input writes 26 symbols 
1 as an output. Set h = niin^ggupp^Q) {b < j}. According to the definition of A, 
for all t €N, /*(c"))h = c^. Hence Condition 2 of Definition [3] is false. Therefore 
if A scattered strictly simulates Ai , 9 has to contain an universal subsequence 
whose support is syndetic. 

On the other hand, assume that 9 contains an universal subsequence a with 
syndetic support. Then, there exists a subsequence a' with a'g < a'^ < . . . < a',^, 
where n is the length of the input of Ai. Build the initial configuration c" 
as described in Construction 3. Since supp(a) has bounded gaps, let p be the 
shortest one and set the radius r — p. The rest of the proof is essentially the 
same as the one given for Theorem [T] with a' playing the role of 6*. □ 

Proposition 2. For any TM Ai that executes in time T{n), consider an ACA 
A given by Construction 3 and an updating sequence 9 containing an universal 
subsequence whose support is syndetic. Then, A scattered strictly simulates A4 
in time O (T(n)2). 



Vi e Z, cf 
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{xi,q,m,2) 
(b, q, m, 2) 
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Proof. Consider A = {A, X,p, 9) where A and A are as in Construction 3. For i G 
N, let Si be as in the proof of Proposition [TJ Let a be the subsequence of 6 given 
in the hypothesis and let p G N be the minimal gap. Consider the subsequence 
a' = S2pS4p . . . S2ip ■ ■ ■ Similarly to the proof of Proposition [U be 
used to encode the simulation of a step of . The total length of the simulation 
is then bounded by 

T(n) T(n) 

^ \s2^p\ + \si,p\ + \se^p\ - ^ 12zp+3 = 6p ( Tiny + (1 + —)T{n) ) - O {T{nf) . 
i=i i=i ^ V J 

□ 



4 Future work 

This paper studies the simulation of Turing Machines by fully asynchronous 
CA. We have shown that computational universality for AC A corresponds to 
the condition that the ACA updating sequence contains any cell infinitely many 
times (universal sequence). We have also exhibited some universal sequences in 
order that the computational cost of the simulation is reasonable (quadratic) 
w.r.t the length of the run of the simulated TM. 

The construction of sequences that allow quadratic simulations might seem 
artificial, classifying our results in the domain of computability but of no use 
in practical simulations. Indeed, for example, consider the program described 
in [T]. It is used to simulate biochemical processes in cells. It is essentially based 
on an ACA (although the authors do not clearly state this) which represents 
proteins (and other chemicals) by particles. The current configuration is updated 
in two steps: diffusion and collision/reaction. In the diffusion step a particle is 
randomly chosen, a direction is randomly chosen and then the particle makes 
a move in this direction. Adopting a more "system" based vision instead of a 
particle based one and passing in ID for simplicity sake, one can represent the 
sequence of activations of the different sites (which may contain a particle or not) 
as 6*0 = 0, 01 = Bq^Xi^Q^ = 6*1 +-^2, ■ ■ ■ ,0t = dt-i+Xt, . . ., where Xi are random 
variables with values in { — 1, +1} identically distributed (with uniform Bernoulli 
distribution for example). Therefore, the current active cell 9t = X]i=o^i 
also a random variable with values in Z and, practically speaking, the updating 
sequence is a random walk (see |10| for more on random walks) . It is well-known 
that ID random walks pass through all sites infinitely many times, therefore 6 is 
an universal updating sequence and hence, by Theorem [T| the system descrived 
in [1] is capable of Turing universal computation. 

For every alphabet A, local rule A and radius r it is possible to consider the 
class of all the ACA such in the form {A, A, r, 6) where 9 is generated by a ID 
random walk. Then it is possible to investigate for any property the probability 
of finding an ACA inside that class that has the given property. In the case of 
simulation of Turing machines it can be proved that for any TM A4 running in 
time T{ri), there exists an ACA with updating sequences generated by random 



walks such that the probabihty of (strictly) simulating Ai in time 3T{n) is 
2-3T(ri) r^YiQ above result says that there exist a class of ACA where a particular 
ACA A chosen uniformly from that class has a low probability of simulating a 
TM faster than deterministic ACA seen in the paper. The authors are currently 
investigating to see if and up to which extent the above results can be improved. 
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